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The temperature dependent evolution of the renormalization effect in optimally-doped Bi2212 
along the nodal direction has been studied via angle-resolved photoemission spectroscopy. Fine 
structure is observed in the real part of the self-energy (ReE), including a subkink and maximum, 
suggesting that electrons couple to a spectrum of bosonic modes, instead of just one mode. Upon 
cooling through the superconducting phase transition, the fine structures of the extracted ReE 
exhibit a two-processes evolution demonstrating an interplay between kink renormalization and 
superconductivity. We show that this two-process evolution can be qualitatively explained by a 
simple Holstein model in which a spectrum of bosonic modes is considered. 

PACS numbers: Valid PACS appear here 



An important issue concerning the microscopic mech- 
anism of high-T c superconductivity is whether Cooper 
pairs are bound by collective modes, analogous to the 
role of phonons in conventional BCS superconductors. In 
conventional superconductors, the electron mass is renor- 
malized by phonons, causing fine structures in the single 
particle spectral function as observed in tunneling spec- 
tra (e.g. Ref. V). In the high-T c superconducting copper 
oxides, an electronic renormalization effect in the form of 
a kink in the band dispersion along the diagonal of the 
Brillouin zone (nodal direction), has been discovered by 
angle-resolved photoemission spectroscopy* 2 * 3 ., suggesting 
that the electrons couple to some bosonic mode(s). Al- 
though the identity of the bosonic mode, lattice^ 5 -^ or 
magnetic ^ 9 ' 10 : 11 , has been the focus of great interest in 
connection to pairing, aspects of this renormalization ef- 
fect still remain unclear. In particular, the lack of a 
superconductivity-induced shift of the energy of the kink 
remains an important puzzle regardless of whether this 
bosonic mode is a phonon or a magnetic mode. 

According to the theoretical calculations using an 
Eliashberg formalism, a coupling to a sharp bosonic mode 
in a d-wave superconductor with a maximum gap Ao 
yields a dispersion kink along the nodal direction at an 
energy of f2 + A q 12 i 13 ' 14 . Above T c , the kink position is 
expected to shift to a lower binding energy due to the 
closing of the superconducting gap; thus, a different kink 
position than that at temperatures below T c is expected. 
This, however, has not been observed; in fact, the "ap- 
parent" kink position in the nodal dispersion remains at 
approximately the same energy (~ 70 meV) in both the 
normal and superconducting state a 8 ' 9 ' 11 . This lack of 
energy shift is a long standing puzzle 15 , which seems to 
contradict theoretical predictions, casting doubts on the 
interpretation of the dispersion kink as a renormalization 
effect due to electron-boson coupling 16 . 



Here, we focus on the band dispersion along (0,0)- 
(7r,7r) direction (nodal direction) and report tempera- 
ture dependent measurements on optimally-doped Bi- 
2212. Compared to other cuprates with lower T c , the 
large superconducting gap (Ao ~ 40 meV) of optimally- 
doped Bi-2212 provides a better opportunity to resolve 
the superconductivity-induced phenomena in the nodal 
electron renormalization. A reference experiment on 



(a), , (b) 




k-k F (n/a) e - E F (eV) 

FIG. 1: (Color online) (a) The MDC-derived band disper- 
sions of optimally-doped Bi2212 at temperature above and 
below Tc. The apparent dispersion kink and the onset of the 
"high energy anomaly" (HEA), are approximately indicated 
by the longer and shorter arrows, respectively. The shaded 
area highlights the difference of band dispersion between these 
two temperatures. The dashed line is the "bare band" for 
extracting the ReE as described in the text. Inset shows the 
MDC-derived band dispersion of non-superconducting Bi2201 
crystals measured at similar temperatures, (b) ReE extracted 
using different choices of the bare band as plotted and labeled 
in the inset. The shaded area indicates the energy of the sub- 
kinks discussed in Fig. [5] 
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heavily overdoped non-superconducting Bi2201 has also 
been performed as a comparison. 

High quality single crystals of optimally doped 
Bi 2 Sr 2 Ca . 92 Y . 08 Cu 2 O 8+5 (B12212 OP, T c = 96K) 17 
were selected to perform the temperature dependence 
experiments. For the reference experiment, sin- 
gle crystals of heavily overdoped non-superconducting 
Pbo.3 8 Bii.7 4 Sri.88Cu0 6 +i (Bi2201 OD, T c < 5K) were 
also measured in the same temperature range. The main 
results of Bi2212 were collected by using He I light (21.2 
eV) from a monochromated and modified Gammadata 
He Lamp with a Scienta-2002 analyzer. The data on 
Bi2201 were collected by using 19 eV photons at Stan- 
ford Synchrotron Radiation Laboratory (SSRL) beamline 
5-4 with Scienta-200 analyzer. The energy resolution is 
~10 meV and angular resolution ~0.35°. The sample is 
cleaved and measured in ultra high vacuum (< 4 x 1CP 11 
Torr.) to maintain a clean surface. The same temper- 
ature dependence experiments on Bi2212 were repeated 
at SSRL, which reproduces the results obtained from the 
He lamp. 

The band dispersion along the nodal direction is ex- 
tracted in the usual manner by fitting the momentum 
distribution curve (MDC) to a Lorentizan function 8 -. The 
band dispersions from above and below T c are shown in 
Fig. [lja). As can be seen, the apparent kink position in 
the dispersion at both temperatures is located at about 
70 meV presenting the aforementioned puzzle . Upon 
close examination, there is a difference between the dis- 
persions above and below T c around the apparent kink 
position (shaded area in Fig. [T](a)). Although this dif- 
ference is subtle, it is significantly larger than that in 
non-superconducting Bi2201 (inset), where the difference 
is due to thermal broadening. This observation implies 
that the dispersion difference seen in the case of opti- 
mally doped Bi2212 is a superconductivity-induced phe- 
nomena, rather than a trivial thermal broadening effect. 
We note that these observations are consistent with the 
previously published data 8 -^^ but the puzzle of missing 
energy shift of the dispersion kink across the supercon- 
ducting phase transition has not yet been addressed. 

This subtle difference could be further analyzed by ex- 
tracting the real part of the self-energy ReE, obtained 
by subtracting the MDC-derived dispersion from an as- 
sumed bare band. A rigorous determination of the bare 
band for this low energy renormalization effect (the kink 
structures) is generally difficult, and further complicated 
by the recently discovered high energy anomaly HEA 18 . 
Nonetheless, to zeroth order, the low energy renormaliza- 
tion effect could be assumed to manifests itself on an ef- 
fective bare band which accounts for all renormalization 
effects from the HEA and other possible contributions. 
We assume this effective bare band can be approximated 
by a linear dispersion connected from the Fermi surface 
crossing to 200 meV below Ep on the dispersion before 
the onset of the HEA, as illustrated by the dashed line 
in Fig. [TJa). We also argue that the temperature de- 
pendence of this bare band is small and negligible in the 
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FIG. 2: (Color online) ReE in the normal state and super- 
conducting state for: (a) optimally-doped Bi2212, (b) Non- 
superconducting heavily overdoped Bi2201. Insets plot ReE 
with their maximum normalized to illustrate the difference of 
the shapes. The arrows indicate the energy positions of the 
kink and subkink. The errors bars were estimated from the 
99.7% confidence interval of the fitted MDC peak positions, 
(c) The evolution of ReE across the superconducting phase 
transition. 



temperature range of interest, since no significant tem- 
perature dependence of the band dispersion is seen at 
energies greater than 150 meV below Ep (Fig. HJa)). 

The extracted ReE of optimally-doped Bi2212 above 
and below T c are plotted in Fig. [U(a). At both temper- 
atures, the maximum of ReE, essentially the apparent 
kink position in the band dispersions, is located at ap- 
proximately the same energy, ^70 meV. This coincidence 
of the maximum positions verifies previous observations 
of the absent shift of the apparent kink position. Impor- 
tantly, an additional feature can now be visualized in the 
extracted ReE from the data taken above T c : a shoulder 
or "subkink" at about ~ 35 ± 4 meV where the slope of 
ReE changes abruptly. We stress that the observation of 
a subkink in ReE above T c in optimally doped Bi2212 
system is a robust experimental feature; it is not only 
well reproducible in different sets of our data, but also 
has been seen in the published data, where this subkink 
feature was not discusse d 11 ! 19 . We also note that while 
the magnitude of the overall self-energy is dependent on 
the choice of the bare band, the position of this subkink 
is relatively immune to this choice as illustrated in Fig. 
QJb). Notably, the fine structures in ReE have also been 
observed in other families of cuprates in the normal state, 
which has been interpreted as the electrons coupled to a 
spectrum of phonons&£. More interestingly, as further 
shown in Fig. [5] (a), comparison of ReE at different tem- 
peratures reveals a clear change of its shape across the 
superconducting phase transition; at a temperature well 
below T c , only a maximum with no "subkink" in ReE 
(dark markers) can be identified. This temperature de- 
pendent shape change is further emphasized in a plot of 
the normalized ReE (inset). 

Further support could be drawn from a reference 
experiment on heavily overdoped non-superconducting 
Bi2201 crystals, whose ReE at two temperatures are plot- 
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FIG. 3: (Color online) (a-b) The self-energy parameter Z{uj) 
of 2 equally strong Einstein modes in the case of (a) nor- 
mal state and (b) a d— wave superconductor with Ao =37 
meV. (c) Extracted ReE by applying MDC-analysis to the 
spectral functions (simulated ARPES data) calculated with 
3 weighted Einstein modes in a tight-bind band structure at 
different temperatures from normal state (104K) to supercon- 
ducting state (88K with A = 20 meV, and 10 K with A = 
37 meV) . The arrows indicate the energy positions of the kink 
and subkink. (d) Band-dispersion along the nodal direction 
extracted from the spectral functions calculated by the same 
condition described in (c). The shaded area highlights the dif- 
ference of band dispersion between these two temperatures. 
Inset plots the MDC-dispersions at 104 K and 10 K extracted 
for the case of a non-superconductor. 



ted in Fig. [3fb). Compared to the data of optimally- 
doped Bi2212 (Fig. Ufa)), there are two major differ- 
ences: (1) no detectable shape change of ReE when 
the temperature varies from 103 K to 10 K; the only 
change here is a slight difference of the magnitude of 
the maximum due the thermal broadening effect, (see 
also the normalized ReE plotted in the inset). (2) No 
prominent subkink structure is seen. These two differ- 
ences strongly suggest that the thermal broadening effect 
can not explain the shape of the ReE and its temper- 
ature dependent evolution observed in optimally-doped 
Bi2212; otherwise, a similar behavior of ReE should have 
been observed in this reference experiment. Therefore, 
what has been observed in optimally-doped Bi2212 is a 
superconductivity-induced phenomena, and the subkink 
in the ReE in the normal state is a generic feature of 
optimally-doped Bi2212. The absence of a subkink in 
heavily overdoped Bi2201 is related to a material and 
doping dependent coupling between the electrons and the 
spectrum of bosonic modes, which has been discussed 
elsewhere^. 

The evolution of extracted ReE from above T c to be- 
low T c in optimally-doped Bi2212 is shown in Fig. [DJc), 
where the data slightly below T c at 88K provides im- 
portant information not available in earlier experiments. 
Two process are involved. First, both the maximum 
(shorter arrow) and subkink (longer arrow), shift toward 



higher binding energy due to the opening of a super- 
conducting gap. This energy shift is expected in the 
Eliashberg formalism lending further support for inter- 
preting fine structures (subkink and maximum) in ReE 
as coupling to a spectrum of bosonic modes. Second, we 
also observed a non-trivial redistribution of the relative 
weight between the lower and higher energy renormaliza- 
tion features across the superconducting transition. The 
lower energy feature (longer arrow) gradually dominates 
in ReE, such that the slope between the two arrows grad- 
ually changes from negative to positive; at a temperature 
well below T c (15 K), the lower energy feature dominates 
and sets the maximum of ReE. The higher energy feature 
is then irresolvable. This gives an impression that there 
is a superconducting-induced change of mode-coupling, 
favoring the lower energy mode at the expense of the 
higher one as the temperature crosses T c . 

To qualitatively understand this spectral weight re- 
distribution of renormalization features in ReE, we con- 
sider a model with a simple bosonic spectrum consisting 
of only two Einstein modes with energies f2 = 36 and 
70 meV, respectively. The mode energies were chosen 
to approximately match the energy of the subkink and 
maximum of our ReE data. The electronic self energy 
ui(l — Z(ui)) is plotted in Figs. [3] (a) and (b) for the case 
of the normal state (i.e. "no gap") and a d-w&ve super- 
conductor with Ao = 37 meV, respectivel y 13 ' 20 . Here we 
take the temperature T=0 and assume a circular Fermi 
surface. In both figures the coupling strength of each 
mode to the electrons is identical (f2 • A = 52 meV is as- 
sumed). Working in weak coupling, the total self-energy 
is a sum of the contributions of both modes as shown in 
the figures. 

When the superconducting gap is absent, the two 
peaks of the total self-energy, renormalization features 
caused by the two coupled modes, are equally strong. 
When the superconducting gap is present, the lower en- 
ergy renormalization feature appears to be stronger than 
the other, even though the coupling strength of these two 
modes are identical. The enhancement of the lower en- 
ergy renormalization feature is essentially what has been 
observed experimentally (Figs. [U (a) and (c)). Accord- 
ing to our model, this enhancement is primary due to 
the highly asymmetric shape of u>(l ~ Z(ui)) near the 
singularity for each mode in the superconducting state; 
it increases monotonically from Ep up to the singular- 
ity fl m ode + Ao, then drops suddenly to a value close to 
zero (Fig. [3Jb)). Thus, when summing self-energies of 
the two modes for the total contribution, the self-energy 
peak induced by the lower energy mode "sits" on the ris- 
ing slope of the self-energy of the higher energy mode; as 
a consequence, it appears to be enhanced; while for the 
higher energy mode, its peak can only sit on the flat tail 
of the self-energy of the lower energy mode, hence, no 
significant enhancement for the higher energy feature in 
the total self-energy. This is different than the "no gap" 
case, where the shape of the self-energy is more or less 
symmetric with respect to the singularity (logarithmic 
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divergence, Fig. [3ja)); thus, the higher energy renormal- 
ization feature also acquires a similar enhancement from 
the tail of the lower renormalization feature. 

The key features demonstrated in Figs. G2a) and (b) re- 
main valid even if the bosonic spectrum consists of more 
than two modes in a generic cuprate band structure at 
finite temperatures. In Fig. Efc-d), we present the simu- 
lated renormalization effect at finite temperatures, using 
the band structure of Ref. ,12, resulting from coupling to 
three modes at energies 36, 50 and 70 meV, correspond- 
ing to energies of typical optical phonons in cuprates. 
The relative strength of each mode coupling has been 
adjusted such that the shape of the ReE in the normal 
state (104 K) looks similar to our 104 K data in Fig. [3] 
(a)2i; the coupling strength then remains unchanged for 
the superconducting state calculations (d-wave gap, at 
88K with A o =20 meV, and 10K with A =37 meV). The 
spectral functions along the nodal direction were calcu- 
lated and then convoluted with a 2D Gaussian response 
function with 10 meV in energy and 0.012 it /a in mo- 
mentum direction to simulate finite instrument resolu- 
tion effect. The MDC analysis used for actual ARPES 
data is applied to this simulated spectral functions to 
extract the renormalized band dispersions. As shown in 
Fig. [31(d), the difference between the MDC-dispersions 
at these two temperatures (shaded area) is significantly 
larger than that caused by thermal broadening (the in- 
set). This reproduces our experimental observations on 
optimally-doped Bi2212 and the reference experiment on 



non-superconducting heavily overdoped Bi2201 (see Fig. 
□ (a)). 

The temperature dependent evolution of the self- 
energy extracted from MDC dispersion is shown in Fig. 
EJc). As expected, the "peaks" in the ReE correspond- 
ing to each individual mode is less resolvable due to the 
superposition with others and the smearing by the reso- 
lution effects, leaving only subkinks in the ReE, reminis- 
cent of the data shown in Fig. [5] and the extracted ReE 
in other cuprates 6,7 . Our simulations also demonstrate 
that the multimode nature can still survive at a tem- 
perature of 104 K in the form of self-energy shape after 
convolving with instrument resolution. Furthermore, the 
two-process evolution of the ReE observed in our data 
(Fig. [2] (c)) is also qualitatively reproduced. We em- 
phasize that the apparent kink (the broad maximum of 
ReE) in the band dispersion at 15 K and 104 K is caused 
by different bosonic modes, not by the same mode. On 
the contrary, in the case of non-superconducting Bi2201, 
the "apparent" kink in the dispersion is caused by the 
same mode, since no temperature dependent energy shift 
of ReE is observed (Fig. [U(b)). These results also warn 
that the "apparent" kink alone does not tell a complete 
story of the renormalization effect in cuprates. 
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